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Abstract 

We derive the set of kinetic equations that control the evolution of gluons in 
the presence of a condensate. We show that the dominant singularities remain 
logarithmic when the scattering involves particles in the condensate. This allows 
us to define a consistent small angle approximation. 


1. Introduction 

In previous papers [T], 0 • the Boltzmann equation including proper treatment 
of Bose statistics, also referred to as the Boltzman-Nordheim equation, was used 
to describe the evolution of a dense system of gluons, such as that produced in 
the early stages of an ultra-relativistic heavy ion collision. The computations 
were done by keeping only elastic scattering in the collision kernel. Also, the 
small scattering angle approximation was used, thereby reducing the Boltz¬ 
mann equation to a Fokker-Planck equation, much easier to solve. A surprising 
conclusion of this analysis was that for a wide variety of initial conditions ap¬ 
propriate for ultra relativistic heavy ion collisions, Bose condensation occurs. 
This is because the density of gluons is too large, with respect to their total 
energy, to be accommodated by a Bose distribution in equilibrium. One might 
object that higher order inelastic scattering terms might alter this conclusion, 
but computations to leading order in inelastic processes suggest that, at least at 
early times, particle number is approximately conserved, and Bose condensation 
is favoredQ (see also discussion in Q and further references therein). 

The analysis to date is incomplete since terms which explicitly depend upon 
the condensate have been ignored in the scattering term of the transport equa¬ 
tion (aside from the inelastic processes just mentioned). Including such con¬ 
tributions is essential to describe the evolution of the system beyond the onset 
of condensation, i.e., with a condensate present. These contributions allow for 
transport of particles into and out of the condensate. It is the purpose of this 
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paper to derive the form of the transport equation in the presence of the con¬ 
densate, in the small scattering angle approximation. In order to do so, we 
need to include the effects of gluon mass and Debve screening in order to regu¬ 
late the infrared and obtain unambiguous resulttlj. The relevant modifications 
of the collision kernel that are due to finite gluon masses are discussed in a 
companion paper [dj • The equations that describe the evolution in presence 
of a condensate have already been obtained in other contexts @, 01 . The 
present paper analyzes the singularities that originate from the scattering of 
zero momentum particles when the interactions are long range, and therefore 
small angles scattering dominate. We show that linear singularities that emerge 
in intermediate stages of the calculation eventually cancel, leaving us with loga¬ 
rithmic singularities only. This allows us to obtain a consistent small scattering 
angle approximation, also referred to as “leading log” approximation. 

The general structure of the transport equation in the presence of a conden¬ 
sate is reviewed in the next section, where the final set of leading log equations, 
which constitute the main result of this paper, are given. In Sect. [3] we present 
the calculation of the individual contributions to the collision kernel, keeping 
in the computation the leading and sub-leading singularities. The cancellation 
of the linear singularities is established in Sect. [H where the consistency of 
the resulting leading log equations is verified by analyzing conservation laws of 
particle number and energy. 


2. The collision term in the presence of a condensate 

We start with the transport equation including only 2 —> 2 elastic scattering, 
i.e. 

If = l f d3p2 d3p3 d3?>4 1 

tJl 2 J {2tt) 3 2E 2 {2tt) 3 2E 3 (2tt) 3 2E 4 2E 4 

x | Mi2^34 I 2 (2^) 4 <5( 4 )(P 1 + P 2 - P 3 - Pa) 

x {////(l + /l)(l + //) - /l/ 2(1 + /3XI + / 4 )} 

(1) 


where the factor 1/2 in front of the integral is a symmetry factor, and ft is a 
shorthand for f p . = f{p^), and similarly, Ei = E p .. Summation over color and 
polarization is performed on the gluons 2,3,4; average over color and polarization 
is performed for gluon 1. We use the small scattering angle approximation for 
the matrix element for gluon-gluon scattering (1 + 2 —>■ 3 + 4), with the dominant 
contributions arising from the t and u exchange channels: 


M 


t s 


_{t — m 2 D ) 2 (u — m 2 D ) 2 


k = 128t r+ IVf 


( 2 ) 


lr The authors of Ref. Q managed, by a suitable choice of the collision kernel, to give a 
meaning to the massless limit. 
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where s, t, u are the standard Mandelstam variables: 

s = (Pi+P 2 ) 2 , t = (P 1 ~P 3 ) 2 , u = (P 1 -P 4 ) 2 , (3) 

and Pi denotes the four-momentum of particle i, Pi = ( E p . ,p/). Momentum and 
energy conservations imply Pl + P 2 = P 3 + P 4 . We have introduced a screening 
mass to_d to regulate the infrared behavior of the one gluon exchange. Note 
that this screening mass is different from the external gluon mass m entering 
the dispersion relation of the colliding gluons, E p = \Jp 2 + m 2 . Both masses 
are supposed to be small compared to the typical momenta of the particles that 
are taking part in the collision. 

Our main goal in this paper is to obtain the form of the kinetic equations that 
govern the evolution of the distribution function when a condensate is already 
formed and participates in the collisions. In this case the distribution function 
is generic.ally of the form 

fip) = n c (27t) 3 <5 (3) (p) + g(p) (4) 

with g(p) referred to as the “regular” part of the distribution. Note however 
that as soon as a condensate is present, g(p) is expected to be singular at p = 0. 
In equilibrium, g{p) would be a Bose distribution with chemical potential equal 
to the mass, i.e., g(p) ~ 2 mT/p 2 at small p. We have verified in a previous 
work that this singular behavior is indeed approached as the system gets close 
to the onset of condensation. We shall assume here that, when the condensate 
is formed, the distribution g(p) keeps a singular component in 1 /p 2 . 

By inserting this two-component distribution back into the transport equa¬ 
tion, we get 

(2tt) 3 S {3) {p^ d t n c + d t g(pi) = C[n c ,g ]. (5) 

Let us examine the collision integral C[n c ,g] on the right hand side of this 
equation. According to Eq. 0 , each distribution /,; in the collision integral 
gives two contributions, one from the condensate, which comes with the delta 
function S < - 3 ' > (p, ; ), and one associated with the regular function g(p i ). 

1. If all four distributions f i j 2 , 3,4 pick up the S 1 - 3 ^ (p,) component, then they 
become identical and the gain and loss terms cancel each other. 

2. The same situation prevails if three distributions pick up the db 3 ) (p 4 ) com¬ 
ponent, since then the fourth must also be at zero momentum due to the 
overall energy momentum conservation. 

3. If two distributions pick up the 6^ 3 \pi) component, then there are two 
situations: (a) one of the factor S^ 3 \p i ) is in the initial state and the 
other one is in the final state —say p 1 = p 3 = 0 in this case energy 
momentum conservation imposes the other two to have equal momenta, 
which results in the cancellation between the gain and loss terms; (b) both 
factors cb 3 ) (pj are in either initial or final states —say p 1 = p 2 = 0 — in 
this case the other two momenta have to vanish too in order to satisfy the 
energy momentum conservation, and we are back to the situation where 
all four momenta vanish, for which the collision term also vanishes. 
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We therefore need to deal only with the terms that contain one or no delta- 
function component. We can write the collision integral as a sum over five 
terms: 

4 

C[n c ,g] = Cb[ff] + y~lCj[n e ,ff]. (6) 

i=1 

The term Cq has no condensate contribution and the remaining four C,’s have 
a condensate contribution, with pi = 0 in C,. 

The next step involves projecting the full transport equation into a set of two 
coupled evolution equations for the condensate and for the regular distribution. 
There is a subtlety here, since, as we shall see, contributions to the condensate 
density can arise from several places. To handle this properly, we define the 
condensate density n c by integrating the distribution function is a small sphere 
of radius po and let po go to zero. That is, 

rPo->° r My, 

n ° = l ” dp J ( 2 # /W ' (7) 

We use this procedure to obtain the equation for the condensate density 

/•Po-*- 0 r -70. 

dtUc = Ci[n c ,g] + J p 2 dp J ^^3 + ^3 + C 4 ] > (8) 

with Ci[n c ,g] = (27r) 3 <5^ 3 )(p 1 ) C\\n c ,g}- While C\ represents an obvious contri¬ 
bution since the delta function is explicit, we shall see that some of the contri¬ 
butions Co, 2 , 3,4 are singular as p x —s- 0 and give also non vanishing contributions 
to the right hand side of this equation. For |p x | > 0, the transport equation 
reads 


dtg(Pi) = C 0 [g} +C 2 [n c ,g} +C 3 [n c ,g] +C 4 [n c ,p]. (9) 

In what follows, we will derive these terms one by one. 

But before we proceed with the calculation of these individual contributions, 
it is useful to present the final result. In the presence of the condensate, the 
kinetic equation gives rise to two coupled equations, which can be written as a 
single equation of the following form (for an isotropic system, with p\ = IpJ) 

^M = -V- l 7 0 (Pi)--V- (10) 

at m 

where the current j7o is a current that depends only on the regular parts of the 
distribution functions and whose expression is identical to what it would be in 
the absence of the condensate (see Eq. (fT5l) below), while J c which multiplies 
n c in Eq. m is given by 

Jc = -Pi Lc [T*g[ + gi (l + 5l )], g[ = ^ (11) 

d-Ei Pi dpi 


4 





and 


A 


TO; 


( 12 ) 



™d 
2 TO 


Finally T* = I a (pi) / If,(pi) is an effective temperature given by the ratio of two 
integrals of the function g(p) (see Eq. (flbl) below). Note that T* is independent 
of pi . The two equations (l53l) and (IMl) can be recovered by separating f(pi) 
into a delta function and a regular contribution, with the procedure indicated in 
Eq. flSJ). These equations are obtained in the leading order of the small scattering 
angle approximation, which consists in isolating the leading singularity of the 
collision kernel as the screening mass tends to zero. As we shall see, there is a 
delicate interplay in the derivation between singularities that arise in the matrix 
element and those that arise form the behavior at small p of the function g{jp). 


3. Determination of the various contributions Ci 

In this section, we derive the expressions of the various contributions Ci to 
the collision integral, and we discuss the elimination of the dominant infrared 
divergences. We consider a system that is uniform and isotropic in momentum 
space, so that the distribution function f(p) depends only on the magnitude 
p = | p\ of the momentum. 

Before we proceed, we note that, except for the contribution Cq which does 
not depend explicitly on the presence of the condensate, for all other contribu¬ 
tions the presence of the condensate has an impact on the form of the function 
g{p) at small momentum. Our analysis relies on the assumption, supported by 
numerical calculations, that, once the condensate appears, g(p) does not deviate 
much from an equilibrium distribution, i.e., g(p) = 1/ {e^ E ^ m ^ T — 1) + Sg. This 
function has a singular behavior ~ 1/e with e = E — m when p —> 0. We shall 
also assume that 5g(p) is regular as p —>• 0, so that for small e, g(p) takes the 
form g(p) ~ T*/e + V*, with T* and V* two constants. 


3.1. The contribution Cq 

The term Cq corresponds to the usual collision term involving no particle in 
the condensate. It is independent of n c and is obtained by replacing all the /j’s 
in the collision integral by their corresponding regular parts gp. 


Co[g} 


j j 

Jp o Jp o 


(2^) 4 <5( 4 )(F 1 +P 2 -P3-P 4 ) 

1 QE\ E 2 E% E 4 . 


M 


ip 2 jp 3 jp 4 

x{ 5354(1 + 5i)(l + 52) - 5i52(l + 53)(1 + 54)}, 


(13) 


where we use the shorthand to denote momentum integration f p = fd 3 p/( 2 n) 3 . 
Note that particle 1 is the particle that one follows in the kinetic equation (i.e. 
we calculate the collision term corresponding to d t fi). 
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The small scattering angle approximation of the above expression has been 
obtained in a companion paper [J]. The result takes the form 


Co[fl] — — V • J'o(pi), (14) 

where the explicit expression of the current is given in . In the very small 
momentum regime p\ <C m, this expression simplifies and reads 

Jolpi^o ->■ -Ko [Ia(pi)g'i +h{pi)gi{l + 9i )], k 0 =36 t xa 2 s C, (15) 


where kq is a constant factor and the two integrals I a {p\) and h(pi) are given 

by 


laipi) 

= El 

/ - 9 ( 1 + 9 ), 

m , 

j p p 

h{pi) 

= El 

[ -(-9\E p )) 


m , 

J P P 


(16) 

(17) 


where g'{E) = dg{E)/dE = ( E/p)d p g. 


3.2. The contribution C\ 

This contribution corresponds to the case where the particle 1 is in the 
condensate (p 3 = 0), and can be kicked out of it by collisions with particle 2. 
It is obtained by replacing, in Eq. CD, fi and 1 + /i by (2?r) 3 <5 (3) (Pi)- We get 

~ _ n Kn c f f 2n5(m +E 2 - E 3 - E 4 ) (E 2 - E 3 )(E 2 +m) 

l[Tlc ' 9i ~ 16 mJ P2 J P3 E 2 E 3 E 4 [(E 3 - to) + m r } 2 

X {5354 - 9293 - 9294 - ff2}, (18) 


where m r = m 2 D /(2m) (see Eq. (fl2]l b Here p 1 = 0, E\ = m, so that s = 
2m(E 2 + to), t = 2m(m — E 3 ), and u = 2m{E 3 — E 2 ). The factor 2 in front 
accounts for the equal contributions from the t and u channels (the collision 
integral is invariant under the exchange of p 3 and p 4 which corresponds to the 
exchange of the u and t channels). In Eq. (1181) . the integration over p 4 has 
been done by using the momentum delta function. At this point, we use the 
remaining energy conservation delta function to carry out the angular part of 
the p 3 -integral. The angular part of p 2 can then be integrated trivially. We get 
(with d Ei = pi&pi/Ei, and p 2 = E 2 - m 2 ) 


k n c I ! 2 j 0(E 2 - E 3 ) (E 2 - E 3 ){E 2 + to) 


Ci1 ’ mI = TSjl l p ‘ dK 


'P 2 1 


’ ( P 2 E 2 ){p 3 E 3 ) [(E 3 - to) + m r f 


x {ff394 - 9293 - 9254 - 92 } 


('E 2 


32tt 3 


dE 2 


dE , 


p\ - ( E 2 + m)[E 3 - to) 
1 Jm [(E 3 -m)+m r } 2 

9394 - 9293 ^ 9294 - 92 f 

J E 4 — yiri-\-E2 — E3 


( 19 ) 
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It is not difficult to see that the Bose-Einstein thermal distribution with g = m 
is still a fixed point, i.e., the combination of statistical factors vanishes for g a 
Bose distribution. The same remark will apply to all the cases to be discussed 
in this section, and will not be repeated. 

In order to obtain the small scattering angle approximation of the expression 
above, we focus on the singular contributions that come from the region £3 > m. 
This is the region from which we expect a singularity in Eq. © when tod —> 0. 
The small scattering angle approximation consists in isolating such singularities, 
and retaining the leading one as the main contribution to the collision kernel. 
As we shall see the leading singularity is a logarithmic one. Other singularities 
emerge, in particular linear singularities, but we shall show later that these 
cancel, leaving us with logarithms only. 

To proceed, we set e = £3 — m and use the notation g e = g(E 3). Then we 
have £4 = to + £2 — £3 = £2 ^ e, and g(E^) ~ (72 — 9 ^ + 5(72 e 2 , where, as 
before, the prime denotes a derivative with respect to the energy. The difficulty 
in the analysis comes from the fact that a small angle scattering involves a 
small momentum transfer, and as a result of such a scattering p% ~ pi = 0. The 
singularity of the matrix element as mj) -> 0, and the presence of a particle in 
the condensate, render therefore the collision integral sensitive to the singular 
behavior at small p of one of the “regular” distribution functions, here g(ps). 
This is the main subtlety in the analysis. Using the notation that we have just 
introduced, we can rewrite Eq. (1191) in the following form: 



( 20 ) 


where A is an ultraviolet cutoff that controls the divergences caused by the 
small e expansion of the distributions (before this expansion, the integrand is 
naturally cut off at energies £ ~ T). By organizing the terms according to the 
various types of e-integrations, we get 


Ci Kij] — ~ 3 — {~H a L g i — Hb Lq + H c L g 2 + Hd L\} 

6 Z 7T m 


( 21 ) 


with the following finite integrals 
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and the following singular integrals 

/■A 


Lgl = 


de 


^9e 


/ o 


l 


L 0 = de 


rjn* 

(e + m r ) 2 m, 

1 1 


+ V * In 


A 

m r 


0 (e + m r ) 2 m r 
A 


L g2 = I de 


e 2 9e 


0 (e + m r ) 2 
A 


T* In 


A 

m r 


L 0 = de 


g (e + m r ) 2 


In 


A 

m r 


(26) 

(27) 

(28) 
(29) 


where we have kept only the dominant terms when tod /A —^ 0. To calculate 
these integrals, we have adopted the general expansion (see the discussion at 
the beginning of this section): g e —>• T*/e + V* + 6(e) and therefore eg e —> 
T* + V*e + o(e 2 ). Note that terms of order 6(e 2 ) or higher are irrelevant. 

The final result for Ci[n c ,g\ can then be written as 


Ci[n c ,g) 


k n c 
32n 3 m 


(- H a T * 



+ (~H a V* + H C T* + H d ) In 



(30) 


3.3. The contribution C 2 

The contribution C 2 corresponds to the case where particle 2 is in the conden¬ 
sate, that is we calculate the collision integral with /2 replaced by ( 2 n ) 3 6 < ' 3 \p 2 ) ■ 
We obtain 


C2[n c ,g] 


k n c 1 f El j pi - {Ex + m){Ej - m) 
167T m p\Ei [{E 4 — m) + m r ] 2 


x 


- 5i53 - 5i54 - 9 \ 


E3—y E\-\- 7 Yl — E4 


(31) 


The dominant contribution in the small scattering angle approximation cor¬ 
responds now to the region E 4 > m. Accordingly, we set e = E 4 — m and 
g e = g(E 4 ). Carrying out the expansion in powers of e as we did earlier, we get 


C 2 [n c ,g] 


k nc 1 f ^ pj ~ {Ex + m)e 
167t m p\Ei J e ^ 0 (e + m r ) 2 

| (32) 


jb(eSeK -5i(l + 5i)] + ' 


(e5e)y +5l5l 


In a similar fashion as we treat C\ above, we can perform the e-integrals and 
simplify the result to be: 


k n c 


167T m 


= 77- -U-K a T*-K b ) — 


m, 


+ (-K a V * + K C T* + K d ) In ( —) 

\m r ) 


( 33 ) 
















with the pi-dependent functions defined as 


K a = P -^l. 

Ei 

K c = — 

E l 


g'{ 


K b = 

Ei± 

pi 


_ pi [51(1 + 51)] 


E ! 


-9 1 


Si 


. Si + TO . 

ffi5i +- 3 —-9i i 1 + 9i) 

Pi 


( 34 ) 


3.^. The contributions 63,4 

Now we move to the cases with either p 3 (C 3 ) or p 4 (C 4 ) picking up the 
delta function component. As we have already noticed, the exchange of p 3 and 
p 4 corresponds to the exchange of the channels u and t , with the symmetry 
C 3 = C 4 and C 3 = C 4 . However, C 3 ^ C 3 , and the same holds for C 4 . This 
is because the contribution to, say, C 3 is dominated by the region p 3 « p 3 , and 
p 3 = 0. As we shall see, a singularity develops then near p 3 = 0. We shall then 
treat separately the two cases C 3 + C 4 and C 3 + C 4 . 


3.4.1. ci + cx 

Let us first compute the C 3 and C 4 contributions following similar proce¬ 
dures as before: 


cl + ct 


k nc 1 f dE pi + (E 2 - m){Ej - E 2 ) 
16tt m piEi J m 2 [{Ei - m) + m r ] 2 


x 


9192 + 5154 + 52ff4 + 34 


E a —y E± +E2 — th, 


(35) 


We proceed to the small angle approximation and introduce the variable e = 
Ei — m = Ei — E% and set g e = g{Ei). We obtain 


cl 


CY 


K n c 1 f dE p\ + (E 2 - m)e 
167T m piEi J m (e + m r ) 2 

x |[(e3 £ )ff2+32(l + 32)]+e 

(36) 


v 92 


( e 9e)'-y + 921 + 92) 


As compared to the previous calculations, the singularity in the limit e —> 0 
(or equivalently p\ —> 0), sits outside the energy integration. To handle this 
singularity, we perform an explicit integration over p 3 in an infinitesimal sphere 
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centered at p x = 0. We get 

k n c 


[ {Cl + CX} = ^ 3 -/ dE\ f 
J P1 ^ o 32 t T 6 mJ E ^ M J A 


dEi I (IE2 

M 


Pl + (E 2 - rn)e 


E^M 

x |[(effe)ff 2 + 52(1 + 32 )] + e 
k n c 


■ + m r y 


,92 


Mf + 92(1 + 92) 


327t 3 m 
k n c 1 


{H a L g i + H b Lq + H e L g 2 + Hf Li} 
1 


327T 3 m (27r 2 ) 


(■ H a T * + H b ) 


m r 


+ {H a V* + H e T* + H f ) In 


_A 

m r 


with 


H P = 


Hf = 


p 2 dE p 


p 2 dE p 


g"(E p ) (E p — to) , 




$'(£ P ) 


(-Ep - to) 

p 2 


g(E p ){l + g(E p )) 


(37) 


(38) 


•(39) 


The result of this integration is finite. This allows us to identify a “hidden” 
singular contribution to the condensate equation, of the form: 


{ci+cn 


(27t) 3 <5 3 (p 1 ) 


k n c 


327T 3 


(H a T* + H b ) 


+ ( H a V * + H e T * + H f ) In (jEj 


(40) 


3.4.2. c% + c\ 

Let us now compute the C 3 and C\ contribution: 


CZ + Cl = 


k n c 1 


dE, 


Pi + {Ei - to) (£4 - Ei) 


167T m piEi J m ‘ [( E 2 — to) + m r ] 2 

x \ - gig2 + gigi + g2gi + gi 


(41) 


E^—tEi+Ez—m 


To proceed to the small angle approximation we introduce the variable e = 
£2 — to = £4 — Ei and set g e = g{E 2 ). We then obtain: 


at + ci = 


k n r 1 


de 


p\ + (Ei - m)e 


167 T m piEi J e _>0 ( e + m r ) 


[{tge)g'i +5i(l +5i)] + e 


( e ffe)^- + g[(l + gi) 


(42) 
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After performing the e-integrals and we get 


ci + c\ 


k n c 
167 t m 


(K a T* + K b ) 



+ (K a V* + K e T* + Kf) In 



(43) 


where the pi-dependent functions are defined as 


Ei 


9i Ei - m , 

t + ^ r 9 \ 


Kf = — 
' Ei 


(i + 9i)g'i 


Ei — m 


pi 


9i(l+9i) 


(44) 


and K a ,K b are given in Eq. (13H) . 


4. Cancellation of Leading Divergences and the Small Angle Equation 

Finally let us put all pieces together. The results that we have obtained in 
the previous section exhibit divergences when one sends the screening mass to 
zero m r -4 0. These divergence are of two types: linear divergence ~ l/m r 
and logarithmic divergence ~ ln[A/m r ]. We shall see that all linear divergences 
cancel out. This leaves only logarithmic singular terms, leading to a consis¬ 
tent small angle approximation for both “ordinary” scattering involving regular 
distribution functions and for scattering to and from the condensate. 


4-1- Cancellation of Leading Divergences 

We first combine the contributions from Ci = (27r) 3 i5 3 (p 1 )C'i and C 3 + Cf, 


as follows: 


Ci + cl + ci 

-»• (2T 

-+ 


= (2tt; 

with 


Hi = 

J p 2 dE p 

h 2 = 

J p 2 dE p 


,3x3, 


32-7T 3 m 


[(- H a T * - H b ) + ( H a T * + H b )\ ( — 


In | A 

m r 


,3x3, 


k n c 1 
^327t 3 m (27r 2 ) 


[HiT* + H 2 ] 


In 


A 

m r 


g"{E p ) + ^ J fg'(E p ) 


= J p 2 dp | d p [p 2 g] | (46) 


9 E 

(1 + 2 g{E p ))g'{E p ) H dLg(E p ) (1 + g{E p )) 


J P 2d P^ d P [p 2 g{E p ){l +g(E p ))\^j . 


(47) 
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We then combine the contributions from C 2 and C 3 + C\, as follows: 


C 2 + Gt + C\ 


— { \(~K a T* - K b ) + (K a T* + K b )] ( — 
1 d 7 t m 1 \m r 

+ [{-K a + K a )V* + (K c + K e )T* + (. K d + K f )\ 


1 


In [ — 

m r 


— — [K\T* 
167 t m 


Ko 


In I A 


with 


K 2 


El 

e 1 

El 

e x 


9 1 


2^1 , 
Pi 91 


~ 2^Pi [pi9l] ) 


(49) 


2 jE 

(1 + 2 gi)g[ - 2 “gi(l + 5i) 

Pi 


1 


= — d Pi [Pidii 1 +ffi)] ■ (50) 
Pi 


dgi.Pi) 

dt 

dn c 

dt 


= {-VJo} 

n r L r 


-L r 


dpi \pi(T* g[ + gi(l + gi))] 


(27 


J p2dp {~^ d p [ p2 ( T *g' + gi 1 + 5 ))] 


1 


rj ^.2 [4*0^0 (-Po)] Iy >0 —^0 


(51) 


JckPx) = —PiL c l^gl+gxil + gx)] . 
With the above, we can then rewrite the equations to be: 


dg(Pi) 

dt 

dn c 

~dt 


= -V- JoiPi) ~-V ■ J c ( Pl ), 

m 

1 KJbM], 1 " e 


27T 2 


Jpo^-0 27T 2 


[Po^c(Po)] 


Po^O ’ 


(52) 

(53) 

(54) 


1 


g(p->0)= g*(p) + 5g(p), g*{p) = e(E _ M)/T , _ 1 


(55) 


(48) 


Leading-log Equations 

Now we can put all pieces together, and obtain the following leading-log 
equations for the co-evolution of the condensate and the regular distribution: 


where L c is given in Eq. m- The new terms associated with the condensate 
can be nicely rewritten after further introducing the following current: 


where it is understood that the first equation is valid for all px > 0 . 

Lastly let us analyze behavior of the Eas. (|53l56[) at extremely small momen¬ 
tum regime. We focus on the evolution after the onset of condensation, i.e. with 
a nonzero condensate n c > 0. As already analyzed previously, the behavior of 
g(p ) at very small p is a local thermal form plus certain corrections, i.e. 
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with T* = I a /Ib and Sg(p ) is supposed to be regular as p — > 0, with Jp(0) finite. 

Let us recall that before onset, n c = 0, and Jo ~ p at small p so that only 
the first equation contributes. Just before the onset the current Jo contains, 
besides the term ~ p at very small p singular contributions in 1 /p and 1/p 2 . 
The singularity in 1/p 2 dominates, but disappears just at onset (when p = m), 
leaving a contribution Jo ~ 1/p. It follows that, after onset, [poJ7o(po)] P o ^ 0 = 
0, and only J c contributes in the equation for n c , i.e., Eq. (15B1) reads 


dn c 

HT 


1 ric 
2ir 2 m 


[pIMpo)] P o ^ 0 . 


(56) 


Since T* = I a /h, the contribution to J c in Eq. (l52l) comes from terms of the 
form Sgigl ~ 1/p 2 (assuming Jpi(0) is finite), so that J c (p — > 0) ~ pq/p 2 , 
which yields indeed a finite contribution to Eq. ([56]) . Note also that while the 
surface term [pq Jo(po)] Po _>o scales as po and hence does not contribute to the 
condensate evolution, it does control the evolution of g[p = 0) in Eq. (l5.il) . since 
the integral of g(p) in a small sphere of radius po also scales as po 

As a consistency check, one can easily verify that the two equations (1551) and 
(1561) satisfy the conservation laws of particle number and energy. Consider first 
particle number conservation: 


dt^J g{p)+n c |=0. (57) 

Because the currents vanish when p oo, the result of integrating the right 
hand side of Eq. (l53l) is 

1 1 n 

^2 [Po^o(po)] po ^ 0 + [p 2 oMpo)\ pq ^ 0 ■ (58) 

The first term vanishes, as we have argued above. The second term cancels 
with the right hand side of Eq. m- For the energy conservation, we need to 
examine the following: 


dt jy E p g(p) + mn c j = J E p [- V • Jo(p)} + ^ j E p [-V • J c 

[Po^c(p 0 )] po ^ 0 • 

The first term on the right hand side may be written as 


E P [—V • Jo(jp)\ = 


f |-V-[E p Jo(p)] + ^- Job)} 
= ^ [pIMpo)\ Pq ^o + J p ^ • Mp) 


= 0, 


(p)] 

(59) 


(60) 


where we have used [poJq(Pq)] P(J ^ 0 = 0^ an d the proof that the second term 
vanishes is identical to that of energy conservation before onset (see paper 1). 
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For the second term of Eq. we have 


He 

m 


E P [-V-J c (p)\ 


Ur 


p 


-f / -V • [EpJ c (p)} + -f- • J c (p) 


to Jp 
m n c 
2tt 2 m 
n c 
2n 2 


[p 2 0 Jc(po)\ 

[pIMpo)] 


E p 

nc 

Po-s-0 m 


Ip E p 


Jc(p) 


Po—>0 ' 


(61) 


where, in the last step, we have used the explicit expression (15121) of the current, 
with T* = I a /Ib and I a , I), defined in Eq. (H6l) . This result is equal and opposite 
to the last term of Eq. (l59l) . which completes the proof of energy conservation. 


5. Conclusion 

The results of our analysis are Eqs. (1531) and (1561) . or equivalently Eq. m- 
We have shown that these equations are of leading order in the small angle ap¬ 
proximation (leading logarithm of the Debye screening mass). These equations 
satisfy energy and momentum conservations, and have a thermal fixed point 
that consists of a finite condensate plus a Bose-Einstein distribution with crit¬ 
ical chemical potential p = to. The behavior of these equations in the vicinity 
of zero momentum was carefully analyzed, exploiting results reported in a com¬ 
panion paper Q , and the main singular behavior of the particle current near the 
origin in momentum space was identified. Numerical analysis of these equations 
will be reported in a later paper. 
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